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The acceleration of charged particles at astrophysical collisionless shock waves is one of the best studied pro- 
cesses for the energization of particles to ultrarelativistic energies, required by multifrequency observations in 
a variety of astrophysical situations. In this paper we discuss some work aimed at describing one of the main 
progresses made in the theory of shock acceleration, namely the introduction of the non-linear backreaction of 
the accelerated particles onto the shocked fluid. The implications for the investigation of the origin of ultra high 
energy cosmic rays will be discussed. 



1. Introduction 

Suprathermal charged particles scattering back 
and forth across the surface of a shock wave gain 
energy. The concept of stochastic energization 
due to randomly moving inhomogeneities was 
first proposed by Fermi |1J. In that original ver- 
sion, the acceleration process is easily shown to 
be efficient only at the second order in the pa- 
rameter V, the average speed of the irregularities 
in the structure of the magnetic field, in units 
of the speed of light. For non-relativistic mo- 
tion, V -C 1, the mechanism is not very attrac- 
tive. The generalization of this idea to the case 
of a shock wave was first proposed in |2I3I4I5| 
and is nicely summarized in several recent re- 
views |ril7l8iyil0| . where the efficiency of the pro- 
cess was found to be now at the first order in 
V. Since these pioneering papers the process of 
particle acceleration at shock waves has been in- 
vestigated in many aspects and is now believed 
to be at work in a variety of astrophysical envi- 
ronments. In fact we do observe shocks every- 
where, from the solar system to the interplane- 
tary medium, from the supernovae environments 
to the formation of the large scale structure of 
the universe. All these are therefore sites of both 
heating of the medium crossing the shock surface 
and generation of suprathermal particles. The 
two phenomena are most likely different aspects 
of the same process, also responsible for the for- 



mation of the collisionless shock itself. One of 
the major developments in the theory of particle 
acceleration at astrophysical shock waves has con- 
sisted of removing the assumption of test particle, 
namely the assumption that the accelerated par- 
ticles could not affect the dynamics of the shocked 
fluid . Two approaches have been proposed to 
treat this intrinsically non-linear problem: the 
two fluid models |llll2ll3ll4H5j and the kinetic 
models |1o117I18I19| , while numerous attempts to 
simulate numerically the process of particle accel- 
eration have also been made |t)W>l2ll^l±-il^4l^| . 
The two fluid models treat the accelerated parti- 
cles as a separate fluid, contributing a pressure 
and energy density which enter the usual conser- 
vation laws at the shock surface. By construction, 
these models do not provide information about 
the spectrum of the accelerated particles, while 
correctly describing the detailed dynamics of the 
fluids involved. The kinetic models on the other 
hand have a potential predictive power in terms 
of both dynamics and spectral shape of the accel- 
erated particles. 

All these considerations hold in principle for 
all shocks but in practice most of the work has 
been done for the case of newtonian shock waves 
(however see for an extension to relativistic 
shocks). Astrophysical studies have shown that 

1 In the presentation at CRIS 2004, particle acceleration at 
relativistic shocks was also discussed at length. That part 
is not included here, due to the limited space available. 
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there are plenty of examples in Nature of flu- 
ids moving at relativistic speeds, and generating 
shock waves. The generalization of the process of 
particle acceleration to the relativistic case rep- 
resents in our opinion the second major develop- 
ment of the theory (Baring, these proceedings). 

In this paper, we will not present a review of 
all the current efforts in the investigation of shock 
acceleration. We will rather concentrate our at- 
tention upon some recent work in the direction of 
accounting for the non-linear backreaction of the 
accelerated particles. 

2. The non-linear backreaction: breaking 
the test particle approximation 

The original theory of particle acceleration was 
based on the assumption that the accelerated par- 
ticles represent a passive fluid, with no dynami- 
cal backreaction on the background plasmas in- 
volved. Within the context of this approximation, 
several independent approaches | llti| give the 
spectrum of the accelerated particles in the form 
of a power law in momentum f(p) = fo{p/po)~ a , 
where the slope a is related in a unique way to 
the Mach number M of the upstream fluid as 
seen in the shock frame, through the expression 
a = AM 2 /(M 2 — 1) (here we asumed that the adi- 
abatic index of the background gas is j g — 5/3). 
This result is easily shown by using the diffusion- 
convection equation in one dimension for a sta- 
tionary situation (namely df/dt = 0): 
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where D is the diffusion coefficient, f(x,p) is the 
distribution function of accelerated particles in 
phase space and Q is the injection function, which 
we will assume to be a Dirac delta function at 
the shock surface x = in the downstream fluid 
(x = + ). The function / is normalized in such a 
way that the total number of accelerated particles 
is given by J£™ x dpAirp 2 f{p). 

As a first step, we integrate eq. ^ around x = 0, 
from x = CP to x = CP, which we denote as points 
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where u± (u 2 ) is the fluid speed immediately up- 
stream (downstream) of the shock and /o is the 
particle distribution function at the shock loca- 
tion. By requiring that the distribution function 
downstream is independent of the spatial coordi- 
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nate (homogeneity), we obtain 

that the boundary condition at the shock can be 
rewritten as 
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We can now perform the integration of eq. |T| 
from x = -co to x — CP (point "1"), in order 
to take into account the boundary condition at 
upstream infinity. Using eq. we obtain 



\p-^-{u 2 - Ui) - ux/o + Qo(p) = 0. 
6 dp 



(4) 



The solution of this equation for fo has the 
form of a power law /o oc p~ a with slope a = 
3ui/(ui — tia) = 3r/(r — 1), where we intro- 
duced the compression factor r = u\jui at the 
shock. For a strong shock r — > 4 and we find 
the well known asymptotic spectrum / — > p~ 4 , 
or N(E) oc E~ 2 in terms of energy (here again 
we assumed that the adiabatic index of the back- 
ground gas is 7 g = 5/3. 

Why should we expect this simple result to 
be affected by the assumption of test particles? 
There are three physical arguments that may 
serve as plausibility arguments to investigate the 
effects of possible backreactions: 1) the spectrum 
E~ 2 is logarithmically divergent in its energy con- 
tent, so that even choosing a maximum momen- 
tum, it is possible that the energy density in the 
form of accelerated particles becomes comparable 
with the kinetic pressure, making the assumption 
of test particles untenable; 2) if the non thermal 
pressure becomes appreciable, the effective adia- 
batic index can get closer to 4/3 rather than 5/3, 
making the shock more compressive and the spec- 
trum of accelerated particles even more divergent; 
3) more divergent spectra imply larger fluxes of 
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escaping particles at the maximum momentum, 
which make the shock radiative-like, again im- 
plying a larger compression and flatter spectra. 

All the three issues raised here point toward the 
direction of making the backreaction more severe 
rather than alleviating its effect, therefore a run- 
away reaction seems likely, which drives the shock 
toward a strongly non-linear cosmic ray modified 
configuration (here the term cosmic rays is used 
in a general way to indicate the accelerated par- 
ticles) . 

We can describe the expected effects on the 
basis of the following simple argument: if, as 
is usually the case, the diffusion coefficient in- 
creases with the momentum of the particles, we 
can expect that particles with larger momenta 
will diffuse farther from the shock surface in the 
upstream section of the gas. At large distances 
from the shock, only the high energy particles 
will be present, while lower energy particles will 
populate the regions closer to the shock surface. 
There is some critical distance which corresponds 
to the typical diffusion length of the particles 
with the maximum momentum achievable, p m ax- 
At this distance, the pressure of the cosmic rays 
is basically zero and the fluid is unperturbed. 
On the other hand, moving inward, toward the 
shock, an increasing number of accelerated par- 
ticles is present, and their pressure contributes 
to the local pressure budget by slowing down the 
fluid (in the shock frame). This effect causes the 
fluid speed upstream to be space-dependent, and 
decreasing while approaching the shock surface. 
The region of slow decrease of the fluid velocity 
is usually called the precursor. The shock, which 
may now be substantially weakened by the effect 
of the accelerated particles, is usually called sub- 
shock. 

It is useful to introduce the two quantities 
Rsub = u\/u2 and Rtot — uo/ u 2, which are re- 
spectively the compression factor at the gas sub- 
shock and the total compression factor between 
upstream infinity and downstream. Here uq, u\ 
and U2 are the fluid speeds at upstream infin- 
ity, upstream of the subshock and downstream 
respectively. The two compression factors would 
be equal in the test particle approximation. For 
a modified shock, Rtot can attain values much 



larger than R su b and more in general, much larger 
than 4, which is the maximum value achievable 
for an ordinary strong non-relativistic shock. 

The shape of the particle spectrum is still de- 
termined by some jump in the velocity field, but 
this quantity is now local: at low energies, the 
compression felt by the particles is ~ R S ub, while 
at p ~ Pmax the effective compression is ~ Rtot- 
It follows that, since R su b < Rtot, the spectrum 
at low energies is steeper than that at higher ener- 
gies: the overall spectrum at cosmic ray modified 
shocks is therefore expected to have a concave 
shape. 

In the following we will describe the effects 
of the particle backreaction following the ki- 
netic semi-analytical approach proposed in |181 
and we will use the most general formalism, 
which includes the possible presence of seed pre- 
accelerated particles in the environment in which 
the shock propagates. We repeat here the steps 
illustrated above for the linear case. Integrating 
again eq. n aroun d x = 0, from x = 0~ to x = + , 
we get Eq. after invoking the homogeneity of 
the particle distribution downstream. 

Performing now the integration of eq. [I] from 
x = — co to x = 0~ we obtain 



1 dfo, , 



uifo + uq/oo + Qo(p)+ 
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Here /oo represents the distribution of seed pre- 
accelerated particles possibly present at upstream 
infinity. 

We can now introduce the quantity u p defined 



Ml 



~T I dx-^f(x,p), 
fo J-oo dx 



(6) 



whose physical meaning is instrumental to under- 
stand the nonlinear backreaction of the acceler- 
ated particles. The function u p is the average 
fluid velocity experienced by particles with mo- 
mentum p while diffusing upstream away from 
the shock surface. In other words, the effect of 
the average is that, instead of a constant speed 
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i*i upstream, a particle with momentum p ex- 
periences a spatially variable speed, due to the 
pressure of the accelerated particles that is slow- 
ing down the fluid. Since the diffusion coefficient 
is in general p-dependent, particles with different 
energies feel a different compression coefficient, 
higher at higher energies if, as expected, the dif- 
fusion coefficient is an increasing function of mo- 
mentum. 

With the introduction of u p , eq. (|SJ) becomes 



1 df 0f , , 



1 du p 
3^ dp 



+u f aD +Q (p) = 



The solution of this equation can be written in 
the following implicit form: 



p #3[u /oo(p) + Qo(p)] 
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In the case of monochromatic injection with mo- 
mentum pi n j at the shock surface, we can write 



The structure of the fluid upstream of the shock 
and the corresponding spectrum of accelerated 
particles is determined if the velocity field U(p) = 
Up/uo is known. The nonlinearity of the problem 
reflects in the fact that U (p) is in turn a function 
of /o as it is clear from the definition of u p . In 
order to solve the problem we need to write the 
equations for the thermodynamics of the system 
including the gas, the rcaccelerated cosmic rays, 
the cosmic rays accelerated from the thermal pool 
and the shock itself. 
0, (7) The velocity, density and thermodynamic prop- 
erties of the fluid can be determined by the mass 
and momentum conservation equations, with the 
inclusion of the pressure of the accelerated par- 
ticles and of the preexisting cosmic rays. We 
write these equations between a point far up- 
stream (x = — oo), where the fluid velocity is uq 
and the density is po = mn gas ,o, and the point 
where the fluid upstream velocity is u p (density 
p p ). The index p denotes quantities measured at 
the point where the fluid velocity is u p , namely at 
the point x p that can be reached only by particles 
with momentum > p. 

The mass conservation implies: 



Qo(P) = — j 2 °{P-Pinj), 
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(9) 



in j 



where n gas ^ is the gas density immediately up- 
stream (x — CP) and rj parametrizes the fraction 
of the particles crossing the shock which are going 
to take part in the acceleration process. 

In terms of R S ub and Rtot, introduced above, 
the density immediately upstream is n gas- i = 
(po/m p )Rtot/Rsub- We can introduce the dimen- 
sionless quantity U(jp) — u p /uq so that 
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(10) 



PoUq = PpUp. 

Conservation of momentum reads: 



(11) 



P UI + Pgfi + P C Rfi = P p U 2 p + Pg^p + PCR, P , (12) 

where P g o and P g , p are the gas pressures at the 
points x = — 00 and x = x p respectively, and 
PcR.p is the pressure in accelerated particles at 
the point x p (we used the symbol CR to mean 
cosmic rays, in the sense of accelerated particles). 
The mass flow in the form of accelerated particles 
has reasonably been neglected. 

Our basic assumption, similar to that used in 
|27j , is that the diffusion is p-dependent and more 
specifically that the diffusion coefficient D(p) is 
an increasing function of p. Therefore the typ- 
ical distance that a particle with momentum p 
moves away from the shock is approximately 
Ax ~ D(p)/u p , larger for high energy particles 
than for lower energy particles 2 . As a conse- 

2 For the cases of interest, D(p) increases with p faster 
than Up does, therefore Ax is a monotonically increasing 
function of p. 
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quence, at each given point x p only particles with 
momentum larger than p are able to affect appre- 
ciably the fluid. Strictly speaking the validity of 
the assumption depends on how strongly the dif- 
fusion coefficient depends on the momentum p. 

The cosmic ray pressure at x p is the sum of 
two terms: one is the pressure contributed by the 
adiabatic compression of the cosmic rays at up- 
stream infinity, and the second is the pressure of 
the particles accelerated or reaccelerated at the 
shock (Pcr(p)) and able to reach the position x p . 
Since only particles with momentum larger than 
p can reach the point x = x p , we can write 



CR,p 



= P 



CR„0 



ICR 



Pcr(p) 



P 



CR 



\u p 



4-7T 
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— / dpp 3 v(p)f (p),(13) 



where v(p) is the velocity of particles with mo- 
mentum p, Pmax is the maximum momentum 
achievable in the specific situation under investi- 
gation, and jcr is the adiabatic index for the ac- 
celerated particles. In Eq. EUthe first term repre- 
sents the adiabatic compression of the pressure of 
the seed particles advected from upstream infin- 
ity, while the second term represents the pressure 
in the freshly accelerated particles at the position 

Xp. 

In the following we use jcr = 4/3 (see JH] 
for a detailed discussion of the reasons for this 
choice). 

The pressure of cosmic rays at upstream infin- 
ity is simply 
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PCR.O = -5- / dpp i v(p)f 00 (p), 



(14) 



where p m in is some minimum momentum in the 
spectrum of seed particles. 

From eq. I|12(l we can see that there is a maxi- 
mum distance, corresponding to the propagation 
of particles with momentum p max such that at 
larger distances the fluid is unaffected by the ac- 
celerated particles and u p = uq. 

We will show later that for strongly modified 
shocks the integral in eq. (|13[1 is dominated by the 
region p ~ p m ax- This improves even more the va- 
lidity of our approximation Pcr, p = Pcr{> p)- 



This also suggests that different choices for the 
diffusion coefficient D(p) may affect the value of 
Pmax, but at fixed p m ax the spectra of the accel- 
erated particles should not change in a significant 
way. 

Assuming an adiabatic compression of the gas 
in the upstream region, we can write 
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where we used mass conservation, eq. The 
gas pressure far upstream is P 9 q = PqUq/ (7 9 Mq ), 
where 7 g is the ratio of specific heats for the gas 
(7 g = 5/3 for an ideal gas) and M is the Mach 
number of the fluid far upstream. 

We introduce now a parameter (;cr that quan- 
tifies the relative weight of the cosmic ray pres- 
sure at upstream infinity compared with the pres- 
sure of the gas at the same location, £cr = 
PcRfil Pgfi- Using this parameter and the defi- 
nition of the function U (p) , the equation for mo- 
mentum conservation becomes 



dU 
dp 



j _ TOR i£R_ U -( 1GR+1) _ I Lr -^ a + i) 



1 dP CR 
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= 0. 
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Using the definition of Pcr and multiplying by p, 
this equation becomes 



dU 
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P 4 v{p)fo{p), 



(17) 



where /o depends on U(p) as written in eq. (|10(l . 
Eq. I|17[) is therefore an integral-differential non- 
linear equation for U (p) . The solution of this 
equation also provides the spectrum of the ac- 
celerated particles. 

The last missing piece is the connection be- 
tween R su b and i?t (, the two compression factors 
appearing in eq. The compression factor at 
the gas shock around x — can be written in 
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terms of the Mach number Mi of the gas imme- 
diately upstream through the well known expres- 
sion 



R 



sub 



fa + mi 

{ lg -l)M* + 2- 



(18) 



On the other hand, if the upstream gas evolution 
is adiabatic, then the Mach number at x — 0~ 
can be written in terms of the Mach number of 
the fluid at upstream infinity M as 



Ml = Ml [ — 



7 9 + l 



= Ml 



7 B +1 



,uq J \ Rtot 

so that from the expression for R su b we obtain 
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Now that an expression between R su b and 
Rtot has been found, eq. l(T7|l basically is an 
equation for i? stl t>, with the boundary condition 
that U(jp max ) = 1. Finding the value of R su b 
(and the corresponding value for Rtot) such that 
U(pmax) — 1 also provides the whole function 
U(p) and, through eq. (JSJ, the distribution func- 
tion /o (p) for the particles resulting from acceler- 
ation and reacceleration in the nonlinear regime. 
When the backreaction of the accelerated parti- 
cles is small, the test particle solution is recovered. 



2.1. Non-linear spectra and the problem of 
multiple solutions 

In Fig. 1 we show an example of the spectrum 
calculated for parameters which are typical of a 
supernova remnant (solid line), as compared with 
the spectrum estimated according with the simple 
model of [23 (broken line) and the result of a 
numerical simulation (dashed line), also reported 
inJ2El- 

In this calculation no seed particles have been 
assumed to be present in the shock environment. 
The good agreement between the semi-analyical 
approach discussed here and the montecarlo sim- 
ulations proves that the semi-analytical approach 
discussed here is quite effective in describing the 
behaviour of cosmic ray modified shock waves as 
particle accelerators. 




10 10* 
p/(m c) 



Figure 1. Example of concave non-linear spec- 
trum obtained with the calculation of |18I19| and 
^cjjescribed here. The parameters used are listed in 
the figure. 



However the situation is in general more com- 
plex than this: previous approaches to the prob- 
lem of cosmic ray modified shocks had already 
shown the appearance of multiple solutions. This 
was first discussed in I12| in the context of two- 
fluid models and in |16I17| by using a kinetic ap- 
proach. Multiple solutions are found with the 
method proposed here as well. In [FiJ\ it was 
pointed out how the multiple solutions appear 
also in the case of reacceleration of seed particles. 
An example of the phenomenon is illustrated for 
the case of no seed particles in Fig. 2, where we 
plot U(p rnax ), bound to be unity for the physical 
solutions, as a function of the total compression 
factor Rtot ■ Here pi n j , p m ax and the shock Mach 
number are all fixed. The solutions are identified 
by the points of intersections of the curves (ob- 
tained for different values of 77, as indicated) with 
the horizontal line at U(p ma x) = 1. 

One can see that for low values of rj (approxi- 
mately unmodified shock) there is only one inter- 
section at Rtot ~ 4. However, when ry is increased 
the intersections may become three. All the three 
solutions are fully acceptable from the point of 
view of conservation laws. For large values of rj 
the shock is always strongly modified (Rtot 3> 4). 
For these cases, the asymptotic shape of the spec- 



7 



100.0 



1 0.0 r 



0.1 



: P =10 5 mc 




max 


7) = 5 V 




/ ?j=10~ 3 | 








-4 \ 







10 



tot 



Figure 2. Example of the procedure that leads to 
multiple solutions. The different curves are ob- 
tained for the values of the parameter r\ indicated 
in the plot. 




0.00 1 0.0 1 0.1 1 10 1 00 1 1 4 1 s 
p/mc 




0.001 0.01 0.1 1 10 100 1000 10" 10 6 
p/ mc 



0.00 1 0.0 1 0.1 1 10 1 00 1 1 4 1 5 

p/ mc 



trum at large momenta is well described by the 
power law f(p) ~ p- 7 ' 2 (or N(E) ~ E^/ 2 ). 

The comparison between the method described 
above and that of |16I17| has been discussed in 
|29j . In Fig. 3, extracted from [2S|, we illustrate 
the spectra and U(p) for a case in which three 
solutions appear (in both approaches). The case 
corresponds to Mach number M = 150, gas tem- 
perature at upstream infinity T — 10 4 K, injec- 
tion momentum pi n j — mc and maximum 
momentum p max = 10 5 mc. In the calculations 
of |16I17) a specific form for the diffusion coeffi- 
cient as a function of momentum is required. For 
reference we adopted a Bohm diffusion coefficient 
D(p) oc p. In Fig. 3, each panel corresponds to 
one solution. We plot in each panel the spectrum 
f(p) multiplied by p (the linear theory would 
predict f(p)p A = constant). The solid lines show 
the spectra as calculated with the approach of 
|18I19| , while the dashed lines are the correspond- 
ing spectra as obtained using the calculations of 
|lf>ll7| with Bohm diffusion. The agreement be- 
tween the two methods is excellent, despite the 
fact that the approach presented here does not 



Figure 3. The figure contains three panels, each 
one illustrating the spectrum and the function 
U{p) for the three solutions obtained using a spe- 
cific set of parameters PU- The solid lines in- 
dicate the results obtained with the calculation 
described here, while the dashed lines are the re- 
sults obtained with the method of |16I17| . 



require the detailed knowledge of the diffusion co- 
efficient. 

2.2. Are multiple solutions related to the 
injection problem? 

The question arises of whether the appearance 
of multiple solutions is an artifact of our igno- 
rance of the parameter m which defines the effi- 
ciency of the shock in injecting particles from the 
thermal pool. Although this is probably not the 
all story, as confirmed by the fact that multiple 
solutions are present even in the case of reaccel- 
eration of pre-accclcratcd particles (in that case 
n is no longer a free parameter) J^j, it is likely 
that injection plays a crucial role. In order to 
show this, we adopt a simple physical recipe for 
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Figure 4. Example of disappearance of multiple 
solutions with the thermal leakage recipe. 



the injection of particles at the shock. Real shock 
fronts are not one-dimensional sheets but rather 
complex surfaces with a typical thickness that for 
collisionless shocks is expected to be of the or- 
der of the Larmor radius of the thermal protons 
downstream of the shock. One should keep in 
mind that the temperature of the downstream gas 
is also affected by the non-linear modification in- 
duced by the accelerated particles, therefore the 
shock is expected to be thinner when the sub- 
shock is weaker. Our recipe is the following: we 
assume that the particles which are injected at 
the shock are those with momentum p > £pth, 
where we choose £ = 3.25 and pth is the momen- 
tum of the thermal particles in the downstream 
plasma (we assume that the gas distributions are 
Maxwellian) , determined as an output of the non- 
linear calculations from the Rankine- Hugoniot 
relations at the subshock. This approximation is 
sometimes called thermal leakage 30 . In physical 
terms, this makes r/ an output of the calculations 
rather than a free parameter to be decided a pri- 
ori. 



In Fig. 4 we plot U(p max ) calculated as de- 
scribed above, in the case in which r\ is evaluated 
self-consistently from the prescription of thermal 
leakage. The different curves are obtained for 
p max = 5 x 10 10 ,10 7 ,10 5 ,10 3 ,10 2 mc (from top 
to bottom) for a fixed Mach number M = 100. 
One can see that only single intersections with the 
horizontal line U(p max ) = 1 are present, namely 
the multiple solutions disappear if the shock is 
allowed to determine its own level of efficiency 
in particle acceleration. This calculation was re- 
peated for different values of the parameters, but 
the conclusion was confirmed for all cases of phys- 
ical interest [2] • One can also see that large val- 
ues of Pmax typically correspond to more modified 
shocks, and that the compression factor can reach 
large values, far from the test particle prediction. 

3. Discussion and General Remarks 

We discussed some aspects of particle acceler- 
ation in astrophysical collisionless shock waves, 
and showed that even when the fraction of parti- 
cles that participate in the acceleration process 
is relatively small (one in 10 4 of the particles 
crossing the shock surface) a large fraction of 
the incoming energy can be channelled into few 
non-thermal particles. This result, found previ- 
ously by using several different approaches, is of 
the greatest importance for the physics of cos- 
mic rays. Not only the accelerated particles can 
keep a substantial fraction of the energy avail- 
able at the shock, but the spectrum of the ac- 
celerated particles may substantially differ from 
a power law, showing a concavity which appears 
to be the clearest evidence for the appearance of 
cosmic ray modified shocks. Despite the passive 
role that electrons are likely to play in the shock 
dynamics, the spectrum of accelerated electrons 
is expected to be determined by the (cosmic ray 
modified) velocity profile determined by the ac- 
celerated hadrons in the shock vicinity. A con- 
cavity in the spectrum of the radiation generated 
by relativistic electrons appears to be one of the 
possible evidences for shock acceleration in the 
non-linear regime. In the case of supernova rem- 
nants, there are hints that this concavity might 
have been observed |32| . 
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One of the aspects of particle acceleration that 
are more poorly understood is the injection of 
particles from the thermal pool of particles cross- 
ing the shock. This ignorance reflects in the diffi- 
culty of determining the fraction of particles that 
takes part in the acceleration process, and we ar- 
gued that this might be the reason (or one of the 
reasons) why calculations of the non-linear shock 
structure may show the appearance of multiple 
solutions. On the other hand, assuming a simple 
recipe for the injection process (thermal leakage) 
is shown to result in the existence of only one so- 
lution. In other words, if the heating and acceler- 
ation processes are interpreted as two aspects of 
the same physical phenomenon, there seem to be 
no ambiguities in the way the shock is expected 
to behave. In this case, there is no doubt that 
strongly modified shocks are predicted. 

The efficient particle acceleration at strong 
shocks is also expected to result in the reduced 
heating of the downstream plasma, as compared 
with the heating achieved in the absence of accel- 
erated particles. This effect should be visible in 
those cases in which it is possible to measure the 
temperatures of the upstream and downstream 
fluids separately, for instance through the X-ray 
emission of the thermal gases. When the shock 
is strongly modified by the accelerated particles, 
a large fraction of gas heating is due to adiabatic 
compression in the shock precursor, rather than 
to shock heating at the gasous subshock. 

In |19| it was pointed out that if the shock prop- 
agates in a medium which is populated by seed 
pre-accelerated particles, the non-linear modifi- 
cation of the shock can be dominated by such 
seeds rather than by the acceleration of fresh par- 
ticles from the thermal pool. This might be the 
case for shocks associated with supernova rem- 
nants, which move in the interstellar medium 
where the cosmic rays are known to be in rough 
pressure balance with the gas. The spectra of 
re-accelerated particles for modified shocks were 
calculated in and showed the usual concavity 
that is typical of cosmic ray modified shocks. 

There is an additional aspect of particle ac- 
celeration at shock waves that has not been dis- 
cussed so far, namely the generation of a turbu- 
lent magnetic field in the upstream section, due 



to the streaming instability induced by the accel- 
erated particles. The fact that the pressure in the 
form of accelerated particles may reach an appre- 
ciable fraction of the kinetic pressure at upstream 
infinity, pou^, suggests that the magnetic field can 
also be amplified to a turbulent value which may 
widely exceed the background magnetic field, and 
approach the equipartition level. In 33 34] the 
process of amplification has been studied numer- 
ically, and this naive expectation has been con- 
firmed. One should however notice that the non- 
linear effects in particle acceleration, discussed in 
this paper, and in particular the spectral modi- 
fication, are not included self-consistently in the 
calculations of the field amplification in the shock 
vicinity. 

All these issues are relevant for the investiga- 
tions of the origin of ultra- high energy cosmic rays 
in many ways: 1) strongly modified shocks can 
be very efficient accelerators, so that the energy 
requirements for the sources we know might be 
substantially relaxed; 2) the spectra of particles 
accelerated at strongly modified shocks are flat- 
ter than those expected in the linear theory. Flat 
spectra generate a GZK feature which is milder 
than that due to steep spectra, therefore it may 
be a less severe problem to explain possible ex- 
cesses of events at the highest energies; 3) mag- 
netic field amplification in the shock vicinity has 
been invoked in the case of SNR's as a possible 
way to accelerate particles up to the ankle in these 
sources [33135136 . For other classes of sources 
this may imply that it is easier to reach ultra- 
high energies in cases that are currently believed 
to have too low magnetic fields. This last point 
deserved deeper investigation. 
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